Reentrant and Isostructural Transitions in a Cluster-Crystal Former 
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We study the low-temperature behavior of a simple cluster-crystal forming system through simu- 
lation. We find the phase diagram to be hybrid between the Gaussian core model and the penetrable 
sphere model. The system additionally exhibits S-shaped doubly reentrant phase sequences as well 
as critical isostructural transitions between crystals of different average lattice site occupancy. Due 
to the possible annihilation of lattice sites and accompanying clustering, the system moreover shows 
an unusual softening upon compression. 
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Since van der Waals and Kirkwood, we better appreci- 
ate the fundamental role of harshly repulsive interactions 
in the organization of matter [1] . But what happens when 
harshness turns into softness? Core softened potentials 
can exhibit microphase separation [3L reentrant melting, 
and isostructural phase transitions[3j , as found in sys- 
tems as diverse as Cerium metal [J|, star polymers [5|, 
dipolar spheres @ , electron bubbles @, Q , and rotating 
Bose gases Q ■ Even softer coreless repulsive interactions 
are also found in complex systems. Nonlinear fields can 
form particle-like structures governed by soliton-like in- 
and for structures with low fractal 
dendrimers 
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teractions 

dimension, such as polymers 



12] 



mm, 



and microgels [15[ , the centers of mass can be immaterial 
and thus overlap with only a finite free energy penalty. 
As a result, materials governed by such interactions ex- 
hibit unusual phenomenology compared to "simple" mat- 
ter 16j. Soft core models are further used to study the 
difficult glass [l7l Hill and classical ground state determi- 
nation problems [l9l liol ] , highlighting the broad interest 
in the behavior of systems with soft potentials. 

A certain universality permeates the thermodynamic 
assembly of systems governed by bounded, purely re- 
pulsive interactions. Two phenomenological categories 
have been identified from a mean-field analysis. Sys- 
tems with pair interactions whose Fourier components 
are purely positive are expected to show reentrant melt- 
ing, while those with some negative Fourier components 
should cluster and freeze into mu ltiply occupied crys- 



tals (MOC) upon compression |17l l2ll427j j . Crossing the 



divide between the two categories can be realized via 
the tuneable generalized exponential model of index n 
(GEM-n) 23], whose pair potential for particles a dis- 
tance r apart is 



-) = ecxp [— (r/a) r 



(1) 



with e and a setting the units of energy and length, re- 
spectively 28]. Upon compression at low temperatures, 
the n = 2 Gaussian core model (GCM) [29j], which has 
a purely positive Fourier transform, forms cubic crystals 



that eventually remelt 3(1 31 1; while all GEM-n of n > 2 
have some negative Fourier components and are thus pre- 
dicted to form MOC at high temperatures [22]], including 
the n — > oo penetrable sphere model (PSM) limit 21]. 

For strong interactions (or effective low temperatures) , 
which may be the experimentally most relevant regime, 
the phase behavior of MOC-forming systems is not un- 
derstood. Various plausible ordering scenarios upon 
compression are suggested by theory and experiments: 
density functional theory predicts a continuous increase 
in clustering [22|, [23|; the PSM limit presents a se- 
quence of second-order phase transitions between crys- 
tals of increasing occupancy [2l|; and bubble solids al- 



ternate liquid and crystal phases of increasing lattice 
occupancy Q. A cascade of pure first-order isostruc- 
tural solid-solid transitions between crystals of different 
occu pan cy is also plausible, as theoretically predicted by 
Ref. [321 ] . Such isostructural transitions should terminate 
at critical points, because the cluster occupancy of high- 
temperature MOC-formers increases continuously with 
density 23, |24j. A rare example of this type of crit- 



ical point is found in pure Cerium, where a pressure- 
induced electronic promotion underlies the transition be- 
tween two isostructural solids with different lattice spac- 
ing 0, [33l ] . Critical points involving volume collapse have 
also been predicted for a varietyof purely classical, dis- 
continuous model potentials pi. l34{ . But the experimen- 
tal colloidal systems in which they could be observed have 
intrinsic size polydispersity [35| or smooth effective inter- 
actions (l6j that inhibit phase separation. The relatively 
broad lattice spacing of MOC and their smooth interac- 
tion form suggest that crystal formation should be less 
sensitive to these experimental constraints. 

In this letter, we present a computational study of the 
low-temperature behavior of the MOC-forming GEM-4. 
Its phase behavior is shown to be a complex hybrid be- 
tween the GCM and the PSM limits. We furthermore 
find complex reentrant transitions and evidence for a cas- 
cade of isostructural transitions. 

We perform lattice Monte Carlo simulations (3(| of 
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FIG. 1: (Color online) T-p (left) and P — T (right) low- and high-T (left inse ts | 24| ) simulation phase diagrams. The isostructural 
critical point (□) is extrapolated from the law of rectilinear diameters [33] and the T — results (•) come from phonon 
theory [32|]. Left: The coexistence regions (shaded) are delimited by simulation results for liquid (A), bcc (V), and fee (0) 
phases. The free energy per particle for the bcc (V) and fee (0) phases at T = 0.03 and p — 0.85 shows that a fee with 
n,c q = 1.94(1) is the ground state (right inset). Right: The phase boundaries (solid lines) are guides for the eye that are 
consistent with the Gibbs-Duhem slopes (not shown) at the coexistence points (■). The triple points (a) are numbered. The 
right inset enlarges the SDR liquid-bccl-liquid-bccl sequence and the bccl-fccl-bccl crystal reentrance. 



the GEM-4 model, whose high-temperature behavior was 
previously determined [24[, for N = 1000-5000 particles 
at constant N, volume V , and temperature T. The pres- 
sure P is obtained from the virial and the Helmholtz 
free energy F of the different phases is calculated via 
thermodynamic integration. The ideal gas reference [37| 
is complemented, for the body-centered cubic (bcc) and 
face-centered cubic (fee) crystal phases, with a poten- 
tial well centered around each of the N c lattice sites [27| ■ 
This crystal reference, by allowing for the characteristic 
multiple occupation of lattice sites and for particle hop- 
ping between those sites, permits a reversible integration 
path. For a fixed number density p, the average lattice 
site occupancy n c = N/N c at equilibrium is identified for 
every state point by simulating crystals at various fixed 
n c then minimizing the resulting constrained free energy 
F(n c ), i.e., identifying the loci F(n° q ) such that (Fig. [1]) 



dF(n c ) 



dn c 



0, 



p,T,n a 



(2) 



and that 



similarly to the scheme employed in Ref. 
used for GEM-4's high-T phase diagram determination. 
In the latter, n c was tuned until the (unphysical) field 



conjugate to n c had vanished 24], 39], which allows for a 
gradient-based minimization of F, but relies on an addi- 
tional independent calculation of the chemical potential 
li. This approach breaks down at low T, where fj, cannot 
be efficiently resolved by Widom's particle insertion 37 1 . 
Phase coexistence is then determined through common 
tangent construction of the free energy data. A linear 
transformation j3F p/N — f3F p/N — up with inverse tem- 
perature {3 and arbitrary unitless parameter k enhances 



the visibility of the coexistence regime (Fig. [2]). 

We present two projections of the the low-T phase di- 
agram in Fig. [TJ As anticipated from the high-T ex- 
trapolation, the cluster bcc phase vanishes at a triple 
point Tj = 0.078(1). But surprisingly the transition is 
preceded by a S-shaped doubly reentrant (SDR) crystal 
phase sequence, where both a maximum and a minimum 
in the coexistence curve are observed. Below the 
phase diagram is surprisingly rich. GCM-like phase be- 
havior 291431 1 is followed by an unexpectedly complex 
clustering regime. Because the tail of the GEM-4 decays 
faster than any inverse power, the liquid freezes into a 
single-occupancy fee (feel) that reaches vanishingly small 
densities at low T^in agreement with predictions from ge- 
netic algorithms [4fj, |4l| and phonon theory [32| . This 
feel phase gives way to a single-occupancy bcc (beel) 

(2) 

phase at a second triple point T t = 0.031(1). For a 

(2) 

narrow temperature range above T t , a beel wedge be- 
tween the liquid and the feci phase leads to reentrant 
crystallization of beel upon compression. The maximum 
freezing temperature for beel at T = 0.039(1) leads to 
a second SDR sequence similar to that observed in the 
hard core/soft shoulder model 0, G|- By contrast to the 
GCM's reentrant melting behavior [29143 lj . here the liq- 
uid reentrance section of the SDR sequence spreads only 
over a finite density regime 0.59 < p < 0.68 and over a 
much smaller temperature range 0.0385 < T < 0.039. 
The intermediate nature of the GEM-4 suggests that 
this reentrance might become more pronounced as the 
GCM is approached, i.e., n — > 2+, and should disap- 
pear before the PSM limit n — > oo, where reentrance 
is not expected [2lj. The connection between the high 
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FIG. 2: Free energy curves of the stable fee (0) and, for 
comparison, the metastable bec (V) structures at T = 0.01 
with k = 142.7. Both curves show the van der Waals loop 
characteristic of a system whose limited size inhibits phase 
separation. The coexistence densities of the feel (left inset)- 
fcc2 (right inset) isostructrual transition are determined by 
the common tangent construction (dashed line). The equilib- 
rium lattice occupancy (•) plateaus near integer values 
for the thermodynamically stable phases. 



and low-T regimes occurs through a third triple point 
T t (3) = 0.040(1), at which beel vanishes. A prior, coarser 
study of the liquid-crystal transition in this regime missed 
both the presence of feel and of the reentrant melt- 
ing [iij]. It also inaccurately assigned the unusual shape 
of the liquid-crystal coexistence curve to the onset of clus- 
tering, while it is rather caused by the reentrant melting. 

Clustering does indeed influence the phase diagram 
topology, but at densities further away from the liquid 
phase. At low T the nature of clustering is unlike what is 
seen at higher T, where n c changes linearly with p result 
ing in a nearly density independent lattice constant 



Here, n c is quasi-quantized, and at very low T the lattice 
constant a changes discontinuously through isostructural 
transitions between fee lattices of nearly perfect integer 
mean occupancy n c «-» n c + 1 (Figs. [2][3]). The first oc- 
currence of these transitions, feel <->■ fcc2, is partially 
interrupted by the beel phase, down to the fourth triple 
point T t (4) = 0.012(1) (Fig. But at higher densities, 
the fcc2 <-> fcc3 coexistence is fully developed. No other 
liquid or crystal phases are found to interfere and genetic 
algorithm results corroborate that no other crystal sym- 



metry should be stable in this density regime [40|, |41 1 . 
It is at the moment computationally difficult to go be- 
yond fcc3, but both a zero temperature treatment paired 
with phonon theory 32| and a simple mean-field cell 
theory predict a cascade of n c <->■ n c + 1 isostructural 
transitions to carry on ad infinitum, slightly broaden- 
ing the coexistence regime between two integer occupan- 
cies. As argued above, the topology of the phase diagram 



FIG. 3: (Color online) Isothermal (T = 0.03) bulk modu- 
lus B (0), its virial contribution (•) and softening correction 
(shaded), and n£ q (V) [IfJ. The solid vertical lines indicate 
phase boundaries. The lattice constant a changes discontinu- 
ously as fcc2— »fcc3 (inset). 



demands that each uninterrupted isotructural transition 
terminates at a critical point, the first one of which is 
found at T c v ' — 0.049(3). Hopping between lattice sites 
should depress j , c (™ c ' ric+1 ) as Uc increases, and mean-field 
critical universality is generally expected 22j,|43j. The se- 
ries of first-order isostructural transitions contrasts with 
the continuous second-order clustering transitions for the 
PSM predicted by cell theory [2lj. This last observation 
suggests that the behavior of the PSM might be singular, 
but further studies are necessary to clarify the GEM-n 
family phase behavior as n — > oo. 

One of the key material properties of MOC is the pres- 
ence of two distinct microscopic mechanisms for respond- 
ing to compression. Like any other crystal, MOC can 
affincly reduce their lattice constant, but additionally 
they can eliminate lattice sites by increasing the mean 
lattice occupancy. We can decompose the bulk modulus 



(3) 



N,T 



into a constant n c virial contribution _B v ; r and a "soft- 
ening" correction B colr , which map directly to the two 
microscopic mechanisms 2J,|44{. At high T, £? corr can 
be as high as half the virial component 24[ , but in low T 
crystals, the quasi-quantized jumps in n c lead to a signif- 
icantly different mechanical behavior. Away from the co- 
existence regions, where n c is nearly constant, the system 
responds only affinely to isothermal compression and the 
virial contribution to the bulk modulus captures the full 
response of the system, i.e., B CO rr ~ 0, (Fig. |3j) [45|. But 
in the softening regions that precede and follow the phase 
transitions the quantization is imperfect, and B C orr ^ 0. 
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Near the bccl-fcc2 transition, for instance, B corr is nearly 
equal to the virial contribution, which means that the 
system exerts almost no resistance to compression. This 
very rapid change in mechanical properties with compres- 
sion is uncommon, and may lead to novel material behav- 
ior. The different physical natures of the virial and soft- 
ening contributions indeed suggest a separation of time 
scales for their microscopic relaxation, with slow parti- 
cle redistributions contrasted by fast affine deformations. 
Hardening or softening of the material upon compression 
might thus depend on the deformation rate. 

We have presented the intriguing low-temperature 
phase behavior of the MOC-forming GEM-4 through a 
numerical method specially designed for this class of sys- 
tems. The complexity of the phase behavior is par- 
ticularly noteworthy considering the simplicity of the 
model, which is free of competing length scales. Exper- 
imental soft matter realizations of MOC are still lack- 
ing, but large-scale, monomer-resolved simulations of am- 
phiphilic cluster-forming dendrimers are currently un- 
der way (4f|. Importantly, the approach outlined in the 
present work should be directly applicable to phenom- 
ena of reversible cluster formation in other branches of 
physics. Examples are: the structures formed by the soft 
solitons [HI 11 1, the quasi-2D electron bubbles in the 
quantum-Hall regime 0, Q , and the predicted clustering 
of vortex lines in rotating Bose gases @ . 
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